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ABSTRACT 

A prominent chemical reaction in interstellar clouds is the formation of molecular 
hydrogen by recombination, which essentially takes place on dust grain surfaces. Ana- 
lytical approaches to model such a system have hitherto neglected the spatial aspects 
of the problem by employing a simplistic version of the sweeping rate of reactants. We 
show how these aspects can be accounted for by a consistent definition of the sweep- 
ing rate, and calculate it exactly for a spherical grain. Two regimes can be identified: 
Small grains, on which two reactants almost surely meet, and large grains, where this 
is very unlikely. We compare the true sweeping rate to the conventional approximation 
and find a characteristic reduction in both regimes, most pronounced for large grains. 
These effects can be understood heuristically using known results from the analysis 
of two-dimensional random walks. We finally examine the influence of using the true 
sweeping rate in the calculation of the efficiency of hydrogen recombination: For fixed 
temperature, the efficiency can be reduced considerably, and relative to that, small 
grains gain in importance, but the temperature window in which recombination is 
efficient is not changed substantially. 
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1 INTRODUCTION 

The complex chemistry of interstellar clouds originates from 
reactions taking place both in the gas phase and on the sur- 
faces of dust grains. The most important surface reaction 
is the formation of molecu lar hydrogen ^ whic h is highly in- 
efficient in the gas phase iVidali et al.ll200ai . Accordingly, 
there has been much recent interest in developing accurate 
and efficient methods for the modelling of chemical reactions 
takin g place on dust grain surfaces (Hcrbst fc Shematovich 
|2003)- The most widespread approach for treating networks 
involving both surface and gas phase reactions is based on 
chemical rate equations, where the populations of differ- 
ent species are described by continuous concentration vari- 
ables and the rate of production of a species is propor- 
tional to the product of th e concentrations of its constituents 
JPickles fc WiUiamsllToTTll . 

However, rate equations become inappropriate for small 
grains subject to low fluxes of reactants, because the con- 
cept of a continuous concentration is questionable when the 
typical number of reactants on a grain d rops below unity 
JCharnlev et al.l Il997t ICaselh et al]ll998l) . In this regime 
it is necessary to consider the full probability distribution 
P{N) of the number of reactants N on the grain (say, 
the number of H atoms), either within a stochastic Monte 



Carlo simulation JTielens fc Hagenlll982l : IStantcheva et all 
I2OOII) or by solving t he master equation governing the time 
evolution of P(N) JGreen et al.ll200ll: iBiham et al.l IJOOll : 
iBiham fc Lipshtatil2002l: iLipshtat fc Bihamll2004^ ■ 

None of the three modelling approaches mentioned 
above takes explicit account of the process by which reac- 
tants migrate and meet on the grain surface, through ther- 
mal hopping or quantum tunnelling. All such spatial aspects 
of the problem are instead lumped into the sweeping rate A, 
which is usually defined as the inverse of the time required 
for a single atom to visit all S adsorption sites on the g rain. 



tor a single atom to visit all o adsorption sites on tne g r 
and is taken to be of the form (IStantcheva et alj|200ir 



A = a/S 



(1) 



where a denotes the rate of hops between neighboring ad- 
sorption sites. This expression is problematic for two rea- 
sons. First, since a random walker in two dimensions oc- 
casionally revisits a site, i t requires more than S hops to 
explore the entir e surface JMontroll fc Weisall965l : iHuehesI 
I1995I : lKruell200a) . Second, and more fundamentally, the ac- 
tual process of interest - the joint diffusion of two atoms 
which terminates when either one of them desorbs, or the 
two meet and react - is distinct from, and not simply related 
to, the exploration of the grain by a single atom. To the best 
of our knowledge, so far the only approach that overcomes 
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these limitations is the fully microscopic Monte Carlo simu- 
lation of the hopping and reaction of mol e cules on the lattice 
of ad sorption sites ijChane et alJ 1200,4 ICxippen fc HerbstI 
120051) . which cannot be used for the modelling of large reac- 
tion networks over long times. 

The goal of this paper is to provide a consistent defini- 
tion of the true sweeping rate, and to show how it can be 
evaluated exactly for a spherical grain. In Sect.|5|we estab- 
lish a transparent interpretation of the sweeping rate and 
express it in terms of the encounter probability of a pair of 
atoms on the grain. This encounter probability is calculated 
for a spherical grain in Sect. El and two regimes of param- 
eters are identified, for which there are simple approximate 
results. This section also highlights a quantitative relation 
between the encounter probability and the deviation of the 
distribution of the number of reactants on the grain from a 
Poisson distribution. In Sect. |1] the results are used to ob- 
tain the true sweeping rate in terms of physical parameters, 
which is then compared to the conventional approximation. 
We find that the true sweeping rate is reduced considerably, 
and the effect is most pronounced for large grains. Random 
walk considerations are employed to reproduce asymptoti- 
cally exact expressions. We then examine the effect of using 
the true sweeping rate on the recombination efficiency (as a 
function of grain size and of grain temperature) in Sect.|S| 
for a set of physical parameters that has been used in other 
work, and that enables us to compare our results to those of 
microscopic Monte Carlo simulations. A summary and our 
conclusions are presented in Sect.|S| 

We focus the discussion on the specific case of hydrogen 
recombination, but expect the results to be relevant also for 
other, more complex diffusion-mediated reactions. 



2 BASIC RELATIONS AND DEFINITIONS 

We consider a grain subject to a flux F of H atoms per unit 
time, which can desorb at rate W . Once formed, molecular 
hydrogen is assumed to desorb rapidly, or at least not to 
affect the diffusion and reaction of H atoms. The standard 
rate equation governing the mean number (A'^) of H atoms 
on the grain then reads JBiham et alJll998r 



= F - W{N) - 2A{NY 



(2) 



d{N) 
dt 

The sweeping rate A appears in the last term on the right 
hand side, which contains the key approximation inherent in 
the rate equation method, i.e. setting the reaction rate pro- 
portional to the product of the mean numbers of reactants. 
Since each reaction consumes two atoms, the recombination 
rate is thus taken to be 



Rr. 



A{NY 



(3) 



Within the master equation treatment, this is replaced 
by the expression JGreen et al.l 1200 ll : iBiham et al.l I2001I: 
iBiham fc Lipshtatil200^ 



Rn 



A{N{N^1)), 



(4) 



which correctly accounts for the fact that the reaction rate 
is proportional to the number of pairs of atoms, and hence 
vanishes when A = 1. Equations ||3J and I^J ^re seen to 
be identical whenever (A^) — (A)^ = (A), i.e. when the 



variance of the distribution P{N) of the reactant number 
equals its mean. This is a defining property of the Poisson 
distribution, which holds irrespective of the mean of the dis- 
tribution. 

The comparison of ^ and (0J shows that a failure of 
the rate equation approach, in the sense of a large discrep- 
ancy between the two expressions, is not necessarily implied 
when the mean number of H a toms (A) becomes small com- 
pared to unity JBiham et alJ =2005'l. The rate equation can 
be safely applied to any arbitrarily small region of a macro- 
scopic surface, as long as atoms can freely enter and exit 
to ensure that the Poisson statistics of the number of re- 
actants in the region is maintained. Rather, the failure of 
the rate equation ^ on small grains is a consequence of 
the confinement of the reactants on the grain surface, which 
implies that every adsorption site is visited many times and 
the probability for the two reactants to meet is sharply en- 
hanced. As a consequence, the probability to find two atoms 
simultaneously on the grain is strongly reduced compared to 
the Poisson distribution, and -Rmastcr ^ Rratccq- A precise 
quantitative relation between the encounter probability for 
two atoms and the deviation of P{N) from a Poisson distri- 
bution will be given below in Sect. 13.41 

The expression (jlj provides a transparent interpreta- 
tion of the sweeping rate A: Since the number of pairs of 
atoms on the grain is A(A — l)/2, 2A is the rate at which 
pairs of atoms are removed by the reaction. The alternative 
pathway for a pair to disappear is through the desorption of 
one of its constituents, which occurs at rate 2W . We con- 
clude that the probability ponc for a pair of atoms to meet 
and react before one of them desorbs is given by 



2A 



A 



pcnc 



2A + 2W A + W 



(5) 



Conversely, the relation ^ can be used to express the sweep- 
ing rate in terms of penc as 



A 



-L ~~ Pcnc 



(6) 



In the next section we show how pcnc can be computed an- 
alytically in terms of the model parameters. 



3 ENCOUNTER PROBABILITY 

Since the confined geometry plays an important role for the 
sweeping rate, we need to specify the shape of the grain. We 
will treat the simplest case of a sphere, which has the crucial 
advantage that all adsorption sites are equivalent. This is in 
contrast e.g. to the disc geome try employed previously in a 
similar calculation llKrual200.'l) . 

The definition of the encounter probability involves a 
pair of atoms, each of which hops to adjacent adsorption 
sites at rate a and desorbs at rate W. Because of transla- 
tional invariance on the surface of the sphere, this situation 
is equivalent to that of a single particle moving with the 
double hopping rate 2a and desorbing at twice the rate 2W, 
while the other atom remains fixed and present all the time. 
As we are only concerned with the encounter probability (as 
opposed to the time of encounter) , multiplying both rates by 
a common factor of 2 does not change the result [compare 
to JSJ]; hence we revert to the original rates a and W for 
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the moving particle and keep the second particle fixed and 
present throughout. 



3.1 Stationary diffusion problem 

Following the approach of iKrud ll2003r . the next step con- 
sists in passing to a continuum limit in which the occupation 
probability n{x) of the moving atom satisfies a stationary 
diffusion equation. The spherical grain has a radius R, while 
the fixed 'target atom' is modelled as a circular surface area 
of radius r, located at the north pole of the sphere. Since 
we assume the grain to be large enough to contain at least 
a few hundred adsorption sites, r ^ _R. In the original dis- 
crete setting the sphere is tiled by adsorption sites forming 
some regular lattice. The radius r of the target atom then 
corresponds to a distance of 2r between adjacent adsorption 
sites. In the appendix^it is shown for some common lattice 
structures that the number of adsorption sites is given by 



S = g{2R/ry, 



(7) 



where the factor g — 0(1) refiects the lattice geometry. 

The diffusion constant D for a two-dimensional random 
walk on the lattice of adsorption sites is related to the hop- 
ping rate a by _D = a{Sx)'^/4, where Sx = 2r is the length o f 
a single step, and hence D = ar^ JMichelv fc Krua 12004) . 
We thus arrive at the stationary diffusion equation 



D\/n + 



F 



^■kR^ 



Wn^Q 



(8) 



describing the occupation probability n of the moving atom 
in the steady state of impingement, desorption and reaction. 
To utilise the azimuthal symmetry of the problem we go to 
spherical polar coordinates, writing the Laplace operator in 
the form 



1 



d 



R^sinOde 



sin 6 



d_ 

do 



(9) 



The reaction removing the pair of atoms is accounted for 
by an absorbing boundary condition at the boundary of the 
fixed atom. 



n{e = r/R) =0. 



3.2 Exact solution 



(10) 



The non-negative solution of ijSl llO^ that is finite at the south 
pole 8 = n is found after transforming to the new variable 
z = cosO, and is given in terms of Legendre functions of the 
first kind Pv{z) by^ 



n{e) 



F 



AtvR'^W 



1- 



P^(-cos6l) 



P.(-cos(r/i?)) 



where the index reads 



i+Vi"(£ 



(11) 



(12) 



^ Mat hematical details can be found in iGradshtevn fc R.vzhilJ 
hOOdt) . sections 8.7 and 8.8, especially 8.706f., 8.823 and 8.840 to 
8.842. Note that, although the index u can become complex, all 
physical expressions are real. 



and we have introduced the diffusion length 
Id = \/D/W. 



(13) 



This is the typical distance an atom on an unbounded sur- 
face would diffuse prior to desorbing. 

The encounter probability for the pair of atoms (already 
averaged over all possible starting points) then consists of 
two parts. The first part is given by the fraction of all im- 
pinging atoms that reach the target by diffusion, that is the 
fraction of the impingement flux F that enters the target as 
a diffusion flux: 



„ / , „N D dn 
Pdiff = 27vRsm{r/R) x "H -^ 



1 = r/R 



F 



Using a well-known identity^ this yields 
sin{r/R) P^ {- cos{r / R)) 



Pdiff 



2(P/£d)2 P4~cos{r/R)) ' 



(14) 



(15) 



where Pj = —(1 — z'^)^'^dPy/d2:. However, a (small) fraction 
of atoms impinges directly on top of the target area, and 
these are responsible for an additional contribution to the 
encounter probability, namely the (purely geometrical) ratio 
« = (1 — cos{r/R))/2 of the target versus the total surface 
area.^ 

The overall encounter probability of a pair of atoms 
therefore is 



Pcnc 



sin{r/R) Pl{~ 



cos(r/P)) 1 



cos(r/P) 



2(P/^D)2P.(-cos(r/P)) 



(16) 



Equation 11611 is the central result of this section. It should 
be emphasised that the encounter probability is independent 
of the impingement flux because it only contains information 
about two atoms that are already assumed to be present on 
the grain. The following section will feature the behaviour of 
the encounter probability in the regimes of physical interest. 



3.3 Limiting behaviour 

We will first derive the approximate behaviour common to 
all relevant regimes (as defined by the ordering of the three 
length scales, r , R and £d) . It was mentioned earlier that r <^ 
R is necessary for the dust grain to host a reasonable number 
of adsorption sites. But clearly, r ^ ^d as well, because 
otherwise the adatom hardly performs any hops during its 
residence time on the grain. In terms of the parameters of 
1161 1 the common feature of all regimes is that r/R <^ 1, thus 
being the appropriate quantity in which to expand first. 

The actual calculation involves some subtleties, and we 
relegate most of this to appendix^ To leading order in r/R, 



2 IGradshtevn fc Rvzhij <2000l) . equation 8.752, 1. 
^ Although this term is of the order of {r/K)^ <^ 1, we found it 
must be accounted for in certain situations. The essential reason 
for this is that, while u <^ Pdift < Pcnc is valid in all regimes 
of interest, we will also be concerned with the complementary 
quantity 1 — pcnc appearing in jSJ. This however can become 
arbitrarily small, and particularly of the order of u or smaller. 
Appendix^will show that inclusion of u even simplifies further 
analysis. 
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one obtains 



-1 

Pcnc 



R 



2 R 

In ( ) - 27- 2'0(!^ + 1) -7rcot(!/7r) 



(17) 

Here tlj{z) = d(lnT)/dz, where r{z) is the F-function, and 
■y = -%I){1) Ri 0.577215665 ... is Euler's constant. To gain 
insight into the behaviour of the encounter probability it is 
useful to examine the two remaining regimes in which r is the 
smallest length scale of the problem. The relation between 
the grain radius R and the diffusion length lu conveys a 
notion of 'small' and 'large' grains. 



3.3.1 Small grains 

Here, the ordering of lengths reads r <^ R <Si io and hence 
\v\ <^ 1. Expanding pcnc about u -^ 0— to first order yields 
Penc ~ 1 + iy[ln{2R/rf - 1], and with ly ^ -{R/luf we 
have^ 



1-4^ 



ln(^ 
r 



-1 



(18) 



In this regime of low desorption rate each atom spends 
enough time on the grain to explore all adsorption sites 
many times. For low fluxes, the recombination efficiency is 
then limited by the rare event that two atoms are simulta- 
neously present on the grain. In this event however, they al- 
most surely meet. It is for encounter probabilities near unity 
that the failure of the rate equation description is most pro- 
nounced, see Sect. 18. 4I 

3.3.2 Large grains 

Now, r < ^D < -R and v = -1/2 -I- i\ with A « R/lu > 1- 
For this form of u, the cotangent in 1171 1 is purely imaginary 
and cancels the imaginary part of ipiy + 1) (cf. footnote 
|1J|. The behaviour of the F-function implies tp[z) « \nz -\- 
0{z~^) for arguments of large modulus, leading to Tietp{v + 
1) « ln(i?/^D) for the remaining real part of tp. Together we 
obtain 



Pcnc 



r) ln(2fe/r)2 



27 



< 1. 



(19) 



For low fluxes this means that the recombination efficiency 
is low because of fast desorption, which does not allow the 
atoms to spend sufficient ti me on the grain to react. T his is 
the s e cond order regime of iBiham fc Lipshtatl (|2002r and 
iKrud (|200.'l ) . in which the confinement of the atoms to 
the grain surface is not felt, and hence the rate equation 
continues to apply even for small numbers of reactants. In 
both limits, l|18|l and 1119^ . the functional dependence on the 
model parameters is the same as that found earlier for the 



disc geometry iKrual2003r . taking into account that there, 
logarithms of large arguments are generally assumed large as 
well (dominating possible numerical constants of the order 
of unity) . 



3.4 Non-Poissonian statistics 

In this section we derive a simple quantitative relation be- 
tween the encounter probabihty and the deviation of P{N) 
from a Poisson distribution, in the limit where the mean 
number of reactants on the grain is small ((A'^) <^ 1). In this 
regime the distribution P{N) can be truncated at A'^ = 2. 
Hence (A^) ^ P{1) + 2P{2) and (A^^) « P(1) + 4P(2), which 
implies that the recombination rate ((IJ is given by 



Rn 



2AP{2). 



(20) 



Alternatively, Rmmstci can be evaluated f or (A'^) <C 1 using 
the fo llowing simple statistical argument iJKrxie et alJl200ri : 
lKrxiBll2003. ') : An atom arriving freshly on the grain will par- 
ticipate in a reaction if, first, another atom is already present 
[which is true with probability P(l)] and, second, if the two 
encounter each other (which is true with probability pcnc). 
Multiplying this with the total flux of atoms onto the grain, 
we have 



Rn 



FP{l)pcnc ~ WP{lfpc 



(21) 



In the second step we have used that, to leading order for 
(A'^) —^ 0, the mean number of atoms on the grain is (A'') « 
F/W « P(l). Combining ^, (|5TJ and ® we arrive at the 
central relation 



P{2)^^il~Pcnc)P{lf 



(22) 



Since for the Poisson distribution P(2) = P{lf/2 for (A^) -^ 
0, we see that the factor 1— Penc is a quantitative measure for 
the deviation of P{N) from the Poisson distribution, which 
reflects the depletion of the probability of pairs caused by 
the recombination reaction. 

The relation 1)22^ translates into a very simple expres- 
sion for the deviation of the rate equation recombination 
rate ^ from the true rate Q. Indeed, for (TV) -^ the 
ratio of the two becomes 



Rn 



^ratccq 



2P{2)/P{1Y^1-Pcnc, 



(23) 



which highlights once more the importance of the confined 
grain geometry in the breakdown of the rate equation de- 
scription. Using 12411 . the ratio 12311 can be directly evaluated 
in terms of the physical parameters of the system. 



■* This approximation does not invalidate the reality of ponc for 
values of u of the form I12i . 

^ We will basically treat logarithms of large quantities as being 
of the order of unity throughout, thence retaining coequal nu- 
merical constants in sums, although the logarithms are usually 
significantly larger. While this might slightly blur the fundamen- 
tal functional relation, it is numerically adequate in many realistic 
situations. 



4 SWEEPING RATE 

We now translate our results back into the original language 
of the discrete picture via the relations {2R/r)^ — S/g, 
£o/r = ^{D/r^)/W = vWW^, and accordingly R/£o = 
^ySW/{4:ga). Rewriting the exact expression p6^ for the 
encounter probability in these quantities is trivial, but only 
useful for the plots to come. 
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4.1 Sweeping rate for small and large grains 

The regime of small grains is now characterized by SW/a ^ 
1, with an encounter probability 

SW 
Pcnc ^ 1 - ^ IHS/g) - 1] , (24) 

and the sweeping rate Q becomes 
W _ -iga/S 



A 



1 - Pcnc ^n{S/g) - 1 



(25) 



For large grains, given by W/a ^ 1 ^ SW/a, we have 
iga 1 



pcnc ' 



SW ln(4a/VK) - 27 ' 
leading to a sweeping rate of 



(26) 



AxWp, 



Aga/S 



(27) 



ln(4a/W^) - 27 ■ 

Compared to the conventional approximation Q , the sweep- 
ing rate is seen to be reduced by a logarithmic factor in both 
limits (the subtracted numerical constants in the preceding 
expressions will be ignored for the remaining discussion, cf 
footnote |5J|. 

In both cases, the argument of the logarithm is of the 
order of the number of all sites that the second (moving) 
atom has visited during its residence time. This behaviour 
is to be expected from random walk considerations for this 
moving atom (we still employ the scenario of the first H atom 
fixed and present throughout, as introduced at the beginning 
of Sect.|21l: It is a known result for two-dimensional random 
walks that the mean number of distinct sites the atom has 
visited after time t (and thus, after at hops on average) is 
(asymptotically) given by 



iVdis 



Tvat 



Cln{Bat) ' 



(28) 



with constants C and B depending o n the type of lattice 
jMontroll fc WeisJl961lHugheJl995^ . Note that this state- 
ment remains true if B is changed to any positive value, and 
accordingly, we will focus on the pre-factor depending on C. 
This constant has a value of C = 1 for a square, C = \/3/2 
for a triangular, and C = 3-\/3/4 for the hexagonal lattice, 
so that for our purposes, we can identify (cf. appendix El 



C = 7r/(4ff) . 



(29) 



The second ingredient of our reasoning is the effective 
residence time trcs of the moving atom in the two limits. If 
the encounter probability is very small (large grains), the 
atom explores only a small portion of the grain, and the res- 
idence time is (to first order in a/{SW)) limited by desorp- 
tion, tics ~ 1/W. The atom has therefore performed a/W 
hops and has visited A'^dis ~ 4:g{a/W) / ln{Ba/W) distinct 
sites. The encounter probability (which includes an average 
over all spatial initial conditions) is then given by the ra- 
tio Ndis/S, reproducing 12611 . This is the probability for the 
atom to recombine before it desorbs, which in this setting 
can as well be obtained as the recombination rate A times 
the residence time 1/W , thereby yielding (I27II once again. 

If, on the other hand, the encounter probability is near 
unity (small grains), the stay of the atom is almost always 
ended by recombination, and the residence time becomes 
fros ~ 1/A. By the same time, the atom has explored the 




logs 



Figure 1. The encounter probability Ponc as a function of the 
grain size S (on a logarithmic scale to the basis of 10). The three 
graphs correspond (left to right) to W/a = 10~^, 10~^ and 10~*, 
and the grain size discriminating between small and large grains 
is then given by the inverse a/W. Starting at unity and monoton- 
ically decreasing, one can see the gentle roll-off of the encounter 
probability to take place around this critical size, yet of very sim- 
ilar shape in all three cases. 



entire grain, so that the number of distinct sites visited has 
saturated at A'dis(tros) ~ S. Inverting this relation one ob- 
tains (again to first order) ires ~ Sln[SB/{4:g)]/{4:ga), thus 
reproducing the true sweeping rate l|25|l . It is also easy to 
regain 12411 without using the general relation @: The prob- 
ability for the atom to desorb before a reaction occurs can 
be written as 1 — pcnc , but in this limit it is also given by 
the desorption rate W times the residence time 1/A. 



4.2 Comparison with the conventional 
approximation 

We now want to illustrate the impact of the reduction factor 
that using the true sweeping rate implies. We will use the 
(appropriately translated) exact expression 1161 for all plots 
throughout. 

Evaluating the reduction factor of the true sweeping 
rate A compared to the conventional choice a/S by ®, we 
get 

A ^ SW Penc ,gQ. 

a/S a 1— Pcnc 

For the rest of this subsection we will not be concerned with 
the influence of the lattice structure, and hence will set g — 
1. With the translation rules put forth at the beginning of 
Sect. 2] one can convince oneself that the only effect that a 
different lattice factor g has on the encounter probability is a 
rescaling of the grain size S — + S/g. It follows that, treating 
the reduction factor 1301 1 as a function of the grain size, it is 
subject only to a rescaling of the variable and the value of 
the function (or a shift on a logarithmic scale, respectively). 

Fig. shows the encounter probability as a function of 
the grain size (on a logarithmic scale to the basis of 10) 
for three representative values of W/a. The ordinate crosses 
the abscissa at 5* « 400 (corresponding to R/r ~ 10) so that 
the region to the right of it is reasonably described by our 
model, cf. Sect. 13. ll 

Fig. 121 presents (with the same conventions and for the 
same values of W/a) plots of the reduction factor A /{a/S) 
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Figure 2. The reduction factor of the true sweeping rate versus 
the conventional approximation, with the same conventions as in 
Fig. Q and for the same values (top to bottom) of W/a. 



as a function of the grain size S. The graphs converge for 
(absolutely) small grains, because upon entering the small 
grain regime, the different parameter W/a does no longer 
affect A/{a/S), cf. I)25|l. With increasing grain size however, 
this reduction decreases to an asymptotic value depending 
on W/a. The approximations given in Sect. l4.l1 for A as well 
as for pcnc are corroborated by these two figures. 

It might come as a bit of a surprise that the reduction 
of the sweeping rate, while of serious account throughout, 
is most pronounced for large, not for small grains. This can, 
however, easily be understood: The key effect for the reduc- 
tion of A as compared to a/S is back-diffusion, the revisit- 
ing of sites by a randomly walking atom. Small grains are 
those which are swept nearly entirely in most cases. This 
means that back-diffusion is rather ineffective in reducing 
the sweeping rate and corresponds to the observation that 
different values of W/a do not change the reduction factor in 
this regime - after all, recombination (not desorption) lim- 
its the residence time, anyway. The larger the grain surface 
becomes, the more does desorption compete with recombi- 
nation as the limiting process. Therefore, there are more and 
more initial conditions for which back-diffusion is effective in 
reducing the sweeping rate. In the limit of large grains, des- 
orption always ends the residence time. As a consequence, 
increasing S no longer increases the total number of sites an 
atom visits, and hence the effect of back-diffusion saturates 
as well, rendering A/ (a/S) constant (as a function of S). 

In Sect. I^^H it was shown that the reduction factor of 
the master equation versus the rate equation recombination 
rate (whenever (A*') <^ 1) exhibits the opposite behaviour 
(it is only relevant for small grains), because this reduction 
is based on a different aspect: The crucial point there was 
the 'non-Poissonianity' of P{N) due to a large encounter 
probability, and ultimately owed to the confinement of the 
reactants to a small surface. 



5 RECOMBINATION EFFICIENCY 

A quantity that is more directly connected to observations 
of diffuse interstellar clouds is the recombination efficiency 
in the production of molecular hydrogen. It should therefore 
prove interesting to examine the effect that the use of the 



reduced, true sweeping rate has on these quantities and to 
compare the results to recent simulations of the process. As 
the underlying analytic model we choose the master equa- 
tion framework that is now generally agreed upon. It is given 
by 



dPjN) 
dt 



F[P{N- l)-P(A)] 



+ W[{N+ 1)P{N + 1)- NP{N)] 
+ A{{N + 2)(A + 1)P{N + 2)- N{N - 1)P{N)] 

(31) 

(with minor modifications for A'^ = 0, D. lGreen et al.lll200lh 
and Biham fc Lipshtat (2002) have independently found the 
analytic stationary solution, yielding 



F Iw/a{2^2F/A) 



(N) = ,, , 

" ^A l^^^_,(2^/2F/A) 

for the moan number of particles on the grain, and 



(32) 



for the recombination efficiency rj = 2Rmastci/F with 
-Rmastcr as iu ||1J . Iv{z) are modified Bessel functions. The 
crucial improvement we introduce is the consistent substi- 
tution of A via JSJ instead of Q. A convenient and phys- 
ically sensible parameter for the impingement flux consists 
in f/W, where / is the (effective) flux per site as defined by 
F — fS. Our full set of parameters now includes W/a, f/W 
and S. 

For easiest comparison with previous work in the field, 
we employ the following standard scenario: We assume des- 
orption and hopping to be thermally activated by the grain 
temperature T, 



W = V exp ( — 



knT 



a = V exp 



kBT 



7^ . (34) 



with a uniform attempt frequency of i^ = lO^'^s"^. In 
an ana lysis of temperatu re-programmed desorption exper- 
iments iKatz et al.l l|l99!Jl found the activation energies for 
an olivine surface to be Eo/ks = 287 K and Ei/ks = 373 K. 
The situation in an interstellar gas cloud with H atoms 
in t he gas phase at a temperature of 100 K is adopted 
from iBiham et al.l (1200 ih and leads to an effective flux of 
/ = 1.8 X 10~®s~^. Altogether the physical parameters are 
thereby given as 



W 



= exp 



)K 



T 



i,= exp (^-47.665). (35) 



5.1 Size dependence 

We will first treat the size dependence of the recombination 
efficiency 77 for given temperature T = 10 K of the grain. One 
can see in Fig. that, first, the true sweeping rate leads 
to a reduced recombination efficiency as compared to the 
conventional approximation in any case. Second, the quick 
decrease in the efficiency as one passes from large grains (in 
terms of a/W) to smaller sizes is shifted to the left, mean- 
ing that now, relative to the overall efficiency, smaller grains 
have gained in importance. For the sake of completeness we 
have further plotted the efficiency as predicted by the rate 
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equation treatment, i.e. JJrateeq = 2_Rrateeq/-F with 7?rateeq 

given by (^J, and (A'^) obtained as tiie stationary solution 
of ||21l . Clearly, using the conventional approximation of the 
sweeping rate then eliminates all dependence of r/ratocq on 
the size of the grain. Choosing the true sweeping rate how- 
ever results in the smaller grains even surpassing the larger 
ones in recombination efficiency, as the size dependence is 
then solely based on that of the reduction factor A/{a/S), 
cf. Fig. 121 This illustrates once more that the replacement 
of the approximate by the true sweeping rate accounts for a 
different aspect of the problem (and is of most concern in a 
different regime) than the transition from the rate equation 
to the master equation treatment, as we already mentioned 
at the end of Sect.0| 

Finally, one can compare our plots with data obtained 
in simulations performed by IChang et alj (|2005|), using a 
continuous-time random walk algorithm. The efficiency we 
obtain for a square lattice is reasonably close to the results of 
their simulations, which backs their assertion that indeed it 
is back-diffusion that significantly undermines the recombi- 
nation efficiency. However, our approach does not reproduce 
the striking difference between a square and a triangular lat- 
tice in the way that the simulations have shown (with the 
triangular lattice giving rise to much more effective recom- 
bination) . 

The origin of this discrepancy is not clear to us at 
present. However, we would like to rule out the possibility 
that it is an inherent shortcoming of the continuum approx- 
imation. While it is true that the transition from a lattice 
random walk to a continuous diffusion picture obliterates the 
fundamental differences of lattice types, many properties of 
a random walk are immune to this process. In particular, we 
have seen in Sect. 14. l1 that the effect of the lattice type that 
is most important to us (the different number of distinct 
visited sites) is accurately represented by the lattice factor 
g in our model. 

As we argued there, the large grain value for the true 
sweeping rate should be given hy A ~ WNdis/S, and the 
residence time is tres ~ 1/W. It follows immediately that 
the asymptotic value of A for a triangular lattice should 
be larger than that for a square lattice by a factor of 
2/-\/3 ~ 1.15, which is the ratio of the respective values of g 
or C~^. Since for large enough grai ns the solution of the mas - 
ter equation yields 77 Ki 2AF/W'^ llBiham fc LipshtadlJOol) . 
this factor should reappear in the efficiencies, which is in 
good agreement with Fig. |3 Judging from the asymptotic 
values for large grains, the marginal difference in recombina- 
tion efficiency between the square and the triangular lattice 
therefore is a reasonable outcome. Genuine lattice random 
walks cannot be expected to yield a significantly larger rel- 
ative difference than our model predicts. 



5.2 Temperature dependence 

The temperature dependence of the recombination efficiency 
is another interesting aspect: Ultimately, from an astrophys- 
ical point of view the puzzle of hydrogen recombination 
in interstellar dust clouds is about t he temperature win - 
dow in which this process is efficient. IChang et al.l 112005^ 
found t hat their simu l ation s gave very similar results com- 
pared to lBiham et alj 11200 jl in the efficiency peak, while for 
higher temperatures the efficiency was somewhat smaller. 



0.15 




0.05 



log 5 



Figure 3. Recombination efficiency »? as a function of the grain 
size S on a 10 K olivine grain in the standard scenario. The thick 
lines show rj{S) as computed with the master equation and the 
true sweeping rate, the continuous line standing for a square, the 
dashed one for a triang ular lattice. The thi ck dotted graph was 
plotted using a/S as in lBiham et alj i200j) . The corresponding 
thin lines show the predictions of the rate equation treatment for 
the respective choice of sweeping rate. 



and again differs for the different lattice types (both effects 
as expected in our model as well, judging from the fixed- 
temperature plot in Fig. I^J. We focus here on the compar- 
ison of the analytical results for the different sweeping rate 
expressions. 

It should be noted first of all that we do not include the 
mechanism of Langmuir-Hinshelwood (LH) rejection. This 
mechanism repels H atoms that impinge on occupied ad- 
sorption sites on the grain surface and is responsible for a 
quick decay of the recombination efficiency for very low tem- 
peratures, where desorption is heavily suppressed and the H 
(and H2) coverage on the grain is large. 

As we neglect this effect and treat a fixed grain size 
(and therefore a fixed impingement flux), the remaining fac- 
tors governing the recombination efficiency are the relation 
between the fixed S and the temperature-dependent W/a, 
and the average number of reactants on the grain, which 
is smaller than, but of the order of F/W . As the tempera- 
ture increases, W/a increases as well and we approach the 
regime of large grains (cf. Sect. 14. it with the resulting small 
encounter probability. Combined with the decreasing F/W 
this prohibits efficient recombination. For lower tempera- 
tures we head towards the regime of small grains, where an 
encounter probability near unity and an increasing number 
of reactants predict a large recombination efficiency. No fi- 
nite temperature completely prohibits diffusion; and due to 
the larger activation energy, desorption is suppressed even 
more. Diffusion will almost always let two atoms meet once 
they are on the surface, and the extremely long time this 
might take at very low temperatures is a hidden feature 
when treating steady-state conditions. 

LH rejection should be negligible for temperatures 
roughly at or large r than the peak efficiency temperature 
JBiham et al.ll200D) . so we can expect our results to validly 
show the breakdown for higher temperatures. A rough esti- 
mate for a temperature below which our results fail can be 
obtained using 13511 . The mean coverage of H atoms on the 
grain is of the order of f/W. The temperature at which this 
coverage becomes of the order of unity surely is beyond the 
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Figure 4. Recombination efficiency as a function of the temper- 
ature T, for a grain of S = 10'^ sites. Again, the continuous line 
shows the master equation result with the true sweeping rate and 
a square lattice, the dashed line represents the result for a trian- 
gular lattice, and the dotted line uses the conventional approxi- 
mation. Note that this plot is horizontally stretched to uncover 
the small differences. The regime of large grains only starts at 
higher temperatures outside the plotted region. 




12 K 



Figure 5. The same plot as in Fig. |1] now for a grain of size 
S = 10^ sites. This size is 'large' for practically all shown tem- 
peratures. 



validity of our assumptions, and this temperature is given 
byrLH«7.8K. 

For a rather small grain of 1000 adsorption sites (Fig.|lJ, 
differences in the various predictions of the recombination ef- 
ficiency are very small, and we increased the resolution by re- 
striction to a very narrow temperature range. The above es- 
timate suggests that the breakdown of rj for T > 8 K should 
also apply to an analogous model incorporating LH rejec- 
tion. On a larger grain with lO" adsorption sites (Fig-EJ one 
can see that recombination ceases to be effective at a lower 
temperature when employing the true sweeping rate rather 
than the conventional approximation (differences between 
the lattice types are still hardly perceptible). That this ef- 
fect only occurs for larger grains is explicable with the help 
of Fig. El which showed that outside the efficiency peak, the 
reduction (of A with respect to a/ S) is most pronounced 
for large grains, since only then back-diffusion becomes ef- 
fective. On the other hand our earlier explanations imply 
that (independent of the fixed grain size) as we increase 
or decrease the temperature, the efficiency has to approach 
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Figure 6. An Arrhonius plot of the recombination efficiency un- 
der the previous conditions. 



the limits zero and unity, respectively. Using the conven- 
tional approximation A = a/ S in yields the quantity 
1/(1 + SW/a) as a substitute for the encounter probability, 
which shares the limiting behaviour of ponc used in the argu- 
ment. Hence in both limits, the efficiency is the same for the 
two sweeping rate concepts. A fairly general suggestion of 
the figures is that grain size hardly seems to matter for the 
upper bound of the temperature range of effective recom- 
bination: A difference of three orders of magnitude of the 
number of adsorption sites results in only a minute change 
of the upper temperature bound. 

Finally, we re-plot the last graph in the Arrhenius fash- 
ion, i.e. as lnJ7(l/r), which is shown in Fig. |S| For rising 
temperature (in the left part of the plot) there exists an 
effective (negative) activation energy, which can be read 
off to be Sact/^B ~ — 440K. Using the large-grain value 
■q ^ 2AF/W^ together with ^ and ^ yields the pre- 
diction {Eq — 2Ei)/kB = — 459K for this energy, in good 
accordance with the estimate. 

For this astrophysical setting, we may sum up the ef- 
fect of employing the true sweeping rate JSJ instead of Q 
as follows: At a given temperature, relative to the decreased 
overall efficiency, smaller grains gain in importance for the 
recombination process, since larger grains are more affected 
by back-diffusion. With a fixed grain size, the correction is 
only pronounced for (absolutely) large grains, but does not 
change the upper temperature bound of efficient recombi- 
nation considerably. This bound rather seems very sensitive 
to a more detailed modelling of the complex s urface struc- 
ture iChang et alJl2005l : IPerets fc Bihamlbooefl . Our results 
thus strengthen the claim that hydrogen formation on dust 
grains crucially depends on the grain surface structure. 



6 CONCLUSIONS 

We could show that it is possible to rigorously define the 
true sweeping rate for diffusion-mediated reactions on a dust 
grain surface, and for a spherical homogeneous grain, it 
could be evaluated exactly. The primary effect of the spa- 
tial structure of the problem consists in a reduction of the 
true sweeping rate compared to the conventional approxima- 
tion, which is due to the back-diffusion of randomly walking 
atoms. This reduction is moderate and dependent on the size 
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of the grain for small grains, whereas for increasing grain size 
it grows stronger before saturating at a value dependent on 
the ratio of the rates of hopping and desorption of reactants. 

The picture of a spherical, homogeneous grain is highly 
idealized, as is the assumption of a single reaction type, but 
since the effects we have found have explanations that re- 
main reasonable for any type of grain and with a very general 
reaction network, we expect the results to remain relevant 
as well. 

Applying the model specifically to the recombination of 
hydrogen atoms on interstellar dust grains, we could com- 
pare the recombination efficiency that our model predicts 
with the results of microscopic Monte Carlo simulations. 
We found that they coincide reasonably well for a square 
lattice of adsorption sites, but that we cannot reproduce, 
and have no explanation for, the substantially increased effi- 
ciency on an (otherwise equal) triangular lattice simulations 
have shown. 

In the wake of the decreased true sweeping rate (com- 
pared to the conventional approximation), the efficiency is 
reduced as well. Least affected by this reduction, small grains 
become more important for the recombination process. For 
a given size of the grain however, the upper temperature 
bound of efficient formation of molecular hydrogen is hardly 
changed. 

For an improved quantitative analysis, the complex sur- 
face structure of interstellar dust grains seems to be of 
paramount importance. It might be a worthwhile enter- 
prise to consider more complicated analytical models which 
can incorporate some of these aspects. A possible next step 
would be the replacement of a single binding energy and dif- 
fusion barrier by energy distributions, accounting for the in- 
homogeneity of the lattice of a dsorption sites on th e grain, as 
has been done in simulations (IChang et al.ll2nn!Th . Another 
interesting appr oach focuses on the eff ect of the porosity of 
the dust grains ijPerets fc Bihamll2nn(tl . 
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APPENDIX A: LATTICE GEOMETRY 

Assuming S to be large, so that curvature effects do not play 
any role for the single atomic 'tile', we can cover a surface 
area 47r_R^ by simply dividing it into such tiles, each of the 
area that a single atom occupies. This tile area depends on 
the distance of (the centers of) nearest neighbor sites, 2r, 
and the lattice type. Elementary trigonometry shows that, 
for a hexagonal lattice (which by common terminology refers 
to a 'honeycomb' lattice formed by the adsorption sites, with 
coordination number 3), it is given by 3y/3r^ , for a square 
lattice it is obviously 4r^, and a triangular lattice (i.e. one 
with coordination number 6) has a tile area of 2\/3r^. So 
we have S = g{2R/r)^ with a factor g = tv/{3\/3) for the 
hexagonal, 7r/4 for the square, and 7r/(2\/3) for the triangu- 
lar lattice, which is smaller than but of the order of unity in 
all three cases. 



APPENDIX B: ASYMPTOTIC EXPANSION 

Asymptotic expansions of the Legendre functions n ear their 
singul arity are given as equations 3.9.1 (13)ff. in lErdelvil 
(Il953l) . Since it was not clear from the outset that the degree 
of these expansions sufficed, we extended the expansion with 
the aid of computer algebra, and en route found a misprint 
(not accounted for in the errata) in equation (15) ibid., which 
erroneously lacks a factor of 2 preceding Euler's constant, cf. 
llE2l f. below. 

We rewrite 11611 as 



Pcnc 



d 



2v{u + 1) dz 



-\nP,{z)] + 



l + z 



(Bl) 



-— cos(r/i?) 

and then use Mathematica to expand the encounter prob- 
ability to first order in (1 -I- z) about 0, keeping u fixed. 
To cast the result into a reasonable form, we manually ap- 
ply various transformations (valid for arbitrary complex v) 
to get rid of all F and most t/j functions, viz. r(l =p i^) = 
±7r/ [sm{uii)V{±u)], i^{±u) = ip{l±u)^:l/v, and ^(iT^^) = 
ipiiiv) ± 7rcot(z/7r). 

Now we identify u = [1 + z)/2, and with the definition 

X ^\nu~'^ -2'y-2ilj{u + l)--KCot{u%) (B2) 

[the square brackets in 117ll 1 we obtain 

—iy(u + l)X —v{i' + l)X'' 

as the expansion of the encounter probability up to linear 
order in u ~ [r/(2_R)]^. The zeroth order term is equivalent 
to equation H17|l and describes a non-trivial behaviour that 
we examine in the main text. The remaining question is why 
the next order is not needed in either of the two regimes 
discussed in Sect. 13.31 As in the main text, we regard the 
lnit~^ term as being of the order of unity (cf. footnote |KJ. 

For large grains, there is no need to include terms of 
0{u): X then is of 0{1) as shown in the main text, while 
~v{v -|- 1) = {R/£d)^ ^ 1. The largest contribution of the 
second term thence has an additional factor oi —v{u + l)u ~ 
[r/{2£-D)]^ <C 1 compared to the first summand. 

The small grain regime requires a bit more effort, since 
1/ and u can generally be of comparable order now. Ex- 
panding about u -^ 0—, we have xf}{l + v) = —7 + 0{i') 
so that X = luM^^ - 1/u + 0{v). Hence -vX = 1 - 
v\nu'^+0{u'^) and {v+l)X = -l/u^\nu-^ ~l + 0[v) = 
-l/v[l-v{\nu'^ -l) + 0{u'^)]. TheOiiP) term thus reads 
1 -I- v{\nu^^ - 1) + 0{v'^). Similarly evaluating the 0{u^) 
terms yields the pre-factor v[l ~u + 0{iy^)]: Without a con- 
tribution of Olv^u^) present, this is but a higher order cor- 
rection to the first summand. One can further see that the 
direct impingement term u included in 11611 is of precisely 
that order, and thence must have gone to cancel a corre- 
sponding term of the expression 1151 . Altogether we have 
shown that 1171 1 of the main text contains all necessary terms 
to consistently perform the subsequent approximations. 
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